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Abstract 


In the UK and other temperate regions, short rotation coppice (SRC) and Miscanthus x giganteus (Miscanthus) are 
two of the leading 'second-generation' bioenergy crops. Grown specifically as a low-carbon (C) fossil fuel 
replacement, calculations of the climate mitigation provided by these bioenergy crops rely on accurate data. 
There are concerns that uncertainty about impacts on soil C stocks of transitions from current agricultural land 
use to these bioenergy crops could lead to either an under- or overestimate of their climate mitigation potential. 
Here, for locations across mainland Great Britain (GB), a paired-site approach and a combination of 30-cm- and 
1-m-deep soil sampling were used to quantify impacts of bioenergy land-use transitions on soil C stocks in 41 
commercial land-use transitions; 12 arable to SRC, 9 grasslands to SRC, 11 arable to Miscanthus and 9 grasslands 
to Miscanthus. Mean soil C stocks were lower under both bioenergy crops than under the grassland controls 
but only significant at 0-30 cm. Mean soil C stocks at 0-30cm were 33.55 + 7.52 Mg C ha! and 
26.83 + 8.08 Mg C ha^! lower under SRC (P = 0.004) and Miscanthus plantations (P = 0.001), respectively. 
Differences between bioenergy crops and arable controls were not significant in either the 30-cm or 1-m soil 
cores and smaller than for transitions from grassland. No correlation was detected between change in soil C 
stock and bioenergy crop age (time since establishment) or soil texture. Change in soil C stock was, however, 
negatively correlated with the soil C stock in the original land use. We suggest, therefore, that selection of sites 
for bioenergy crop establishment with lower soil C stocks, most often under arable land use, is the most likely 
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to result in increased soil C stocks. 
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Introduction 


Tackling climate change is one of the greatest challenges 
facing the world (IPCC, 2014). Along with other renew- 
able energy sources and demand reduction, the use of 
biomass as a low-carbon (C) replacement for fossil fuels 
is seen as an essential part of the move towards a more 
sustainable energy system (Renewable Energy Road 
Map 2007; DECC et al., 2012). Sources of biomass are 
diverse and include waste streams from food, forestry 
and conventional agricultural crops (Rowe et al., 2009; 
DECC et al., 2012). There is, however, increasing interest 
and utilisation of so-called second-generation (2G) 
bioenergy crops, especially in temperate developed 
nations such as Europe and the USA (Davis et al., 2012; 
Don et al., 2012). These 2G bioenergy crops, predomi- 
nantly perennial grass and woody species, are grown 
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specifically to use as a renewable fuel source and are 
characterised by low input requirement and high 
growth rates. These traits result in a low energy require- 
ment per unit of energy produced, limited management 
requirements, potentially higher C savings and reduced 
environmental impacts when compared to conventional 
food crops used for the production of first-generation 
biofuels (Fazio & Monti, 2011; Don et al., 2012; Mohr & 
Raman, 2013; Walter et al., 2014). 

Assessing the C balance of 2G bioenergy crops pre- 
sents a unique challenge as, in contrast to the use of 
conventional agricultural crops or waste streams, 
bioenergy crop production requires a major change in 
land use and management (Rowe et al., 2009; Aylott & 
McDermott, 2012; Mohr & Raman, 2013). Land-use 
change (LUC) is known to be a primary factor affect- 
ing soil C stock (Guo and Gifford, 2002), and whilst 
impacts of harvesting and utilisation of these crops on 
the C balance are relatively well understood, impacts 
on soil C stocks are less well defined (Fazio & Monti, 
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2011; Rowe et al., 2011; Don et al., 2012; Walter et al., 
2014). 

In their meta-analysis, Don et al. (2012) highlighted 
the limited number of studies on the impacts of bioen- 
ergy crops on soil C stocks in temperate regions, and 
the highly variable and sometimes contradictory results 
reported across these. Even within single multi-site 
studies, impacts on soil C stock have been found to be 
variable between sites, with Walter et al. (2014), for 
example, reporting rates of change in soil C stocks 
across 21 SRC plantations in central Europe from —1.3 
to 1.4 Mg C ha! yr! for transitions from arable land 
and —0.6 to 0.1 Mg C ha! yr for transitions from 
grassland. Meanwhile, for Miscanthus transitions from 
arable land, Poeplau & Don (2014) found rates of 
change in soil C stocks within their study ranging from 
—0.17 to 1.54 Mg C ha! yr! and ranges in the litera- 
ture of between —6.85 and 4.51 Mg C ha! yr 

Some of the variations in the observed impact on soil 
C stocks, both between and within studies, have been 
related to differences in climatic conditions, original 
land use, soil types, management or crop genotype 
(Don et al., 2012; Poeplau & Don, 2014; Richter et al., 
2015). These sources of variability can help to improve 
understanding of the mechanisms underlying changes 
in soil C stock, but comparison of studies can also be 
confounded by differences in quantification methods 
(Don et al., 2012; Bárcena et al., 2014). For example, LUC 
to SRC and Miscanthus can result in changes in soil C 
distributions within the soil profile and therefore sam- 
pling depth, which often differs between studies, can 
have a profound effect on the quantified impacts on soil 
C stocks (Poeplau & Don, 2014; Walter et al., 2014). In 
their meta-analysis of impacts on soil C stocks of LUC 
to forestry, Bárcena et al. (2014) also highlighted the fail- 
ure of many studies to adjust for change in soil bulk 
density (BD) that often co-occur with LUC. This results 
in an incorrect assessment of change in soil C stock and 
inflated between-study variability (Bárcena et al., 2014). 
Apart from some notable exceptions (Walter et al., 2014; 
Ferchaud et al., 2015), few temperate bioenergy LUC 
studies have directly addressed the issue of changing 
BD (Don et al., 2012). 

In the context of mainland GB, and for the two domi- 
nant bioenergy crops in the UK, SRC willow and Mis- 
canthus (Aylott & McDermott, 2012), we address these 
issues by providing a methodologically consistent data 
set of the impacts on soil C of land-use transitions to 
these crops, whilst incorporating variability in potential 
regulatory factors such as climate. This study aims both 
to assess within mainland GB the current impacts on 
soil C stocks of LUC to commercial plantations of either 
SRC or Miscanthus, and to provide insights and data on 
regulatory factors that can be incorporated into future 


modelling activities (see Dondini et al., 2015). To meet 
these aims, we undertook the assessment of soil C 
stocks under 20 Miscanthus and 21 SRC commercial 
plantations and their paired controls. Transitions were 
located across mainland GB and were purposefully 
selected to cover a wide range of climatic and soil con- 
ditions, including soil texture, pH, initial soil C stocks, a 
range of bioenergy crop ages and land-use transitions 
from both grassland and arable land uses, thus allowing 
the influence of these factors on changes in soil C stocks 
to be explored. Soil sampling utilised a combination of 
0-30-cm and 0-1-m soil cores and soil C stocks were 
adjusted for changes in bulk density. 


Materials and methods 


Site selection 


A database of potentially suitable commercial SRC and Mis- 
canthus plantations was populated through liaising with bioen- 
ergy companies and individual growers. Data on soil C stocks 
prior to the land-use change were not available for these com- 
mercial sites, thus a paired-site approach was utilised, where 
impacts on soil C stock are assessed through a comparison 
between a target land use and an adjacent paired control repre- 
senting the original land use (Davis & Condron, 2002; Laga- 
niere et al., 2010). The paired-site method assumes no pre- 
existing differences between the control and bioenergy land 
uses that would confound changes in soil C stock (Wellock 
et al., 2011; Hewitt et al., 2012). Bioenergy plantations were 
therefore selected on the basis of the availability of a suitable 
paired control field in addition to the bioenergy crop age (time 
since establishment), geographical location and the type of 
LUC (ie. from arable or from grassland). Selection aimed to 
provide the widest range of bioenergy crop age and geographi- 
cal location, and a balance of transitions from arable and grass- 
land to SRC and Miscanthus (Table 1). Each control and 
bioenergy plantation pair is referred to as a transition. In total, 
41 transitions were assessed at 28 locations across mainland GB 
(Fig. 1). 

The 41 transitions comprised 12 arable to SRC (all willow), 9 
grasslands to SRC (8 willows, 1 poplar), 11 arable to Miscanthus 
and 9 grasslands to Miscanthus transitions (Table 1). Grassland 
was defined here using Defra definitions and includes both 
permanent pasture (>5 years old) and temporary grassland 
(5 years old and under), with the majority of sites being perma- 
nent pasture (Table 1). The lower number of grassland transi- 
tions reflects the greater difficulty experienced in locating 
bioenergy plantations established on former grassland. 


Sampling method 


Surface soil (0-30 cm). The surface soil of the cropped area of 
each bioenergy plantation or control field was sampled using a 
hierarchical design (Keith et al., 2014), developed to capture 
variability across different spatial scales (Conant & Paustian, 
2002; Conant et al., 2003). Five sampling plots per field were 
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Fig. 1 Map of sampling locations. Dark grey = SRC willow, 
light grey = Miscanthus; the data points of different bioenergy 
crops present at the same location are offset. 


randomly selected from intersections of a grid overlaid on a 
map of the cropped area of field. The resolution of the grid 
was adjusted to ensure that there were a minimum of 50 grid 
intersections, with the condition that the resolution of the grid 
could not be <5 m. A 20-m perimeter buffer was also used to 
reduce potential edge effects. Within the five sampling plots, 
the three within-plot soil cores were taken using a split-tube 
soil sampler (Eijkelkamp Agrisearch Equipment BV, Giesbeek, 
The Netherlands) with an inner diameter of 4.8 cm to a depth 
of 30 cm. The first core was taken at the grid intersect, with 
two further cores taken at distances of 1 m and 1.5 m in ran- 
dom compass directions from the intersect. This gave a total of 
15 spatially nested samples per field, accounting for both field- 
scale (between sampling plots) and plot-scale (cores within 
plots) variability. Before each core was taken, litter (L) and fer- 
mentation (Là horizons were collected from a 25 cm x 25 cm 
area centred on the coring location. Soil cores were divided in 
the field into 0-15 cm and 15-30 cm (measuring from the base 
of the core), individually bagged and returned to the labora- 
tory. There was limited compression in some cores and this 
was allocated to the 0-15 cm section under the observation that 
most compression occurred in the upper layer of soil. The 
depth of the hole was always measured to ensure that the accu- 
rate core length was known. 


Deep cores (0—100 cm). One of the five sampling plots was 
randomly selected and three 1-m cores were taken following 
the same spacing as the 30-cm cores, with the exact coring loca- 
tions adjusted to avoid those of the 30-cm cores. Cores were 
taken using a window sampler system with a 4.4 cm cutting 
diameter (Eijkelkamp Agrisearch Equipment BV, Giesbeek, The 
Netherlands), allowing a full 1-m core to be extracted and sub- 
sequently transported in one section. If coring to the full depth 
was not possible, for example when large stones or bedrock 
were encountered, the precise depth of the cored hole was 
recorded. 


Laboratory processing 


Litter samples were dried at 80 °C for 24 h and dry mass of 
woody material (e.g. twigs, branches), leaves and undifferenti- 
ated material was recorded. Litter was assumed to have C con- 
centration based on litter dry mass of 43% and 45% for 
Miscanthus leaves and stem, respectively (Beuch et al., 2000; 
Robertson et al, in preparation), 42% and 49% for willow 
leaves and stems, respectively (Chauvet, 1987; Heller et al., 
2003), 46% for grass litter (Ross et al., 2002) and 41% for cereal 
litter (Aita et al., 1997). 


Short cores (0-30 cm). The fresh mass of the 0-15 cm and 
15-30 cm core sections was recorded and sections were then 
cut lengthways into quarters for separate subsequent analy- 
ses. One quarter was then set aside for processing for soil C 
and bulk density (BD, Table S1), together with the large 
stones and roots (25 mm) hand-sorted from the remaining 
three sections. Another quarter was used to assess soil pH 
(Table $1) and the remaining sections were archived as a fro- 
zen sample (—20 9C). 

For the assessment of soil pH, the fresh samples were bulked 
within each sampling plot but not across depths giving 10 com- 
posite samples per site (five each for the 0-15 cm and 15- 
30 cm depths). The fresh, bulked samples were sieved to 4 mm 
to remove stones and roots. 10 g of bulk soil was then mixed 
well with 25 ml of deionised water and allowed to stand for 
30 min, before the pH of the liquid layer was recorded (Hanna 
pH210 Meter, Hanna Instruments Ltd., Befordshire, UK). 

For BD, texture and soil C assessment, the fresh soil mass 
was recorded and then samples were air-dried at 25 °C for a 
minimum of 10 days. Air-dried samples were reweighed, 
sieved to 2 mm and the mass and volume of stones and 
roots remain on the sieve recorded. A subsample of the 
sieved soil (15-18 g) was oven-dried (105 °C for 12 h) and 
moisture-loss was recorded. The oven-dried subsample of soil 
was grounded in a ball mill (Fritsch Planetary Mill) and a 
100-mg subsample was used for the assessment of C concen- 
tration using an elemental analyser (Leco Truspec CN, Milan, 
Italy). Prior to analysis using the elemental analyser, soil sub- 
samples that were either from sites located on soil types 
known to contain inorganic C or which had pH values > 6.5 
were tested for the presence of inorganic C using acid fumi- 
gation following Harris et al. (2001). All samples from sites 
which tested positive were treated to remove inorganic C fol- 
lowing the same procedure. 
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A subsample of the sieved air-dried soil was also used to 
assess soil texture. As for pH measurement, samples were 
bulked across each field but not across depth, thus giving one 
value per field for each depth (0-15 cm and 15-30 cm). Analy- 
sis of the bulked samples was conducted by Macaulay Scien- 
tific Consulting Ltd. (Aberdeen, Scotland) with proportions of 
sand, silt and clay analysed by laser diffraction (Malvern 
Mastersizer 2000, Malvern Instruments Ltd., Worcestershire, 
UK). Analysis was conducted for both the bioenergy crops and 
the paired controls. 

Bulk density of the whole core was calculated using values 
of moisture-loss from the air and oven-dried subsamples fol- 
lowing methods in the GB Countryside Survey (Emmett et al., 
2008; Reynolds et al., 2013). These calculations accounted for 
the measured mass and volume in the soil cores taken up by 
stones, and so are corrected to represent the fine earth propor- 
tion (Schrumpf et al., 2011). The Countryside Survey conducted 
a pilot study to compare different protocols to estimate BD in 
different soil types and found that the method used in this 
study was consistent with other protocols and within the 
ranges of typical values expected for each of the soil types 
(Emmett et al., 2008). 

The soil C concentration and bulk density data were used to 
derive mass-based values of soil C stock to account for differ- 
ences in bulk density across transitions. A soil C stock was cal- 
culated based on an equivalent soil mass approach (ESM), 
using a reference dry soil mass of 3 Gg ha 1. following the 
method of Gifford & Roderick (2003). 


Deep cores (0—1 m). On return to the laboratory, the 1-m 
cores were divided into three sections: 0-30, 30-50 and 50- 
100 cm. In cases where compression of the core had occurred 
during sampling, the length of the sections was reduced to 
account for the compression; a method also utilised by Wal- 
ter et al. (2014). Depth increments of 0-30 cm, 30-50 cm and 
50-100 cm were selected based on the common use of these 
increments in similar LUC studies (Laganiere et al., 2010; Don 
et al., 2012). 

Each 1-m core section was divided lengthways, one-half, 
and all root and stones (25 mm) were processed for bulk den- 
sity and C content as outlined for the 30-cm surface soil cores. 
The remaining half was retained as a frozen archive. 

Soil C stocks were again calculated based on an equivalent 
soil mass approach (ESM), using a reference dry soil mass of 6 
and 13 Gg ha^! for the 0-50 cm and 0-1 m sections, respec- 
tively, following Gifford & Roderick (2003). 


Treatment of under length core 


In the ESM calculation, the length of the cores is not directly 
used to calculate soil C stocks (a reference mass is used and 
the deepest sections are used only to give C concentration). It 
is still necessary, however, to remove from the data set any 
cores that, due to the present of large stones or bedrock, do 
not reach a depth that provides a representative C concentra- 
tion for the deeper soil layers. Therefore, based on inspection 
of the soil C profiles, cores «22.5 cm and 70 cm in length for 
the 30-cm and 1-m cores, respectively, were removed from 


the data set prior to statistical analysis (see Table S2 for 
details). 


Statistical analysis 


The difference in soil C stock and litter variables between the 
land uses (SRC, Miscanthus, arable and grassland) was tested 
using linear mixed-effect models with the nlme package in the 
R statistical program (Pinheiro et al., 2014). Differences were 
observed in the control fields of the bioenergy crops with a 
higher overall mean soil C in the arable control sites of the 
Miscanthus transitions compared to the SRC transitions. The 
inclusion of site as a random factor was not sufficient to 
account for this underlying bias and, consequently, the SRC 
and Miscanthus transitions were analysed separately. Land use 
was entered as a fixed effect and field nested within site and 
plot nested within field entered as random effects in all mod- 
els to ensure that appropriate comparisons of transition units 
were accounted for within site. The significance of the vari- 
able land use in the model was examined using a likelihood 
ratio test compared to the null model, including only random 
terms. 

The significance of differences between the levels within 
‘land-use’ was tested using Tukeys multiple comparison in the 
glht function in the multcomp package (Hothorn et al., 2008). 
Marginal (R?) and conditional (R2) R? values were calculated 
(Nakagawa & Schielzeth, 2013; Johnson & O'Hara, 2014) using 
the r.squaredGLMM function (Lefcheck, 2014) in the MuMIn 
package (Barton, 2015). Data on soil C for ESM at 0-30 cm and 
0-1 m were log-transformed prior to testing to meet model 
assumptions. Litter data were x * 1 log-transformed due to 
high number of zero values in arable control fields. In all cases, 
means and standard errors given for land-use effects refer to 
model-estimated values, and therefore account for the random 
effect of site. 

Difference in mean soil C stock between the controls and 
their paired bioenergy crops was divided by the age of the 
bioenergy plantation to estimate annual rates of change in soil 
C as Mg C ha ! yr". This procedure standardises differences 
in soil C stocks between the SRC and the Miscanthus control 
fields, allowing SRC and Miscanthus transitions to be combined 
into the same statistical test. Differences in annual rates of 
change between the 4 transitions (arable to SRC, grassland to 
SRC, arable to Miscanthus and grassland to Miscanthus transi- 
tions) were tested using a two-way ANOVA, with fixed factors of 
control land use (grassland or arable) and bioenergy crop (SRC 
and Miscanthus). Site was not included as a random factor as it 
was not found to improve the model fit. 

Linear regression focused on the 0-30 cm depth where 
change was most likely and was used to explore variables influ- 
encing the impacts of transition to bioenergy crops on soil C 
stocks (clay content, soil pH, soil C stocks, bioenergy crop age, 
MAP and MAT). Data on percentage change from control were 
tested, again to standardise differences in soil C stocks between 
SRC and Miscanthus control fields, allowing SRC and Miscanthus 
transitions to be combined into the same statistical test. 

The drivers of soil C changes were identified through 
model selection but the number of data points limited the 
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complexity of candidate models. Therefore, R? (Nakagawa 
& Schielzeth, 2013; Johnson & O'Hara, 2014) and Akaike's 
Information Criterion (AIC) were first used to assess the 
influence of each explanatory variable on the percentage 
change in soil C stock (Table S3). The explanatory variables 
were then added consecutively to the final models in the 
order indicated by greater R7 or lower AIC scores, pro- 
vided the AIC of model continued to decrease. Site was 
included as a random variable in each model and calcula- 
tion was performed in R using the r.squaredGLMM (Lef- 
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check, 2014) and AIC functions in the Lme4 and MuMIn 
package (Barton, 2015; Bates et al., 2015). 

Selection based on both the R2, and the AIC scores resulted 
in the selection of the same model which included the fixed 
factors control soil C stock and the bioenergy crop type and the 
random effect of site (Table S4). The significance of the 
explanatory variables within this model was examined using a 
likelihood ratio test. 

Prior to this analysis, exploration of the soil texture data 
showed that in contrast to the percentage sand and silt, which 
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Fig. 2 Control versus bioenergy crops soil C stocks for the SRC transitions: 0-30 cm (a) 0-100 cm (c) depths, and the Miscanthus tran- 
sitions 0-30 cm (b) 0-100 cm (d) depths; red symbols represent ex-arable transitions, and green symbols represent ex-grassland transi- 


tions. * indicates site 17 vs. 17C. Error bars give standard error. 
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showed a correlation between the bioenergy crops and paired 
control (R? = 0.72 and R? = 0.71, respectively; Fig. S1), the cor- 
relation for the percentage clay content was poor (R? = 0.26; 
Fig. S1). This poor correlation appeared to be related to high 
soil inorganic C in five of the sites (2, 5, 7, 17 and 19), a factor 
known to affect laser assessment of clay content (Kerry et al., 
2009). Removal of these sites resulted in an improvement to an 
R? of 0.62 but did not improve the explanatory power of the 
percentage clay in regard to the percentage change in soil C 
stocks (Table S3). Thus, this subset was not used in any subse- 
quent analysis (Table S3). 


Results 


Soil C stocks 0-30 cm 


Land use was found to affect surface (0-30 cm) soil C 
stock. (Mg C ha ! in both the SRC (%™(3) = 15.30, 
P = 0.001, R2—0.86) and Miscanthus transitions Es 
(3) = 13.71, P = 0.001, R? = 0.92) (Fig. 2a,b). The greatest 
differences in soil C stocks were in the grassland transi- 
tions, with mean soil C stocks 33.55 + 7.52 Mg C ha^! 
and 26.83 + 8.08 Mg C ha! lower under the SRC 
(P = 0.004) and Miscanthus plantations (P = 0.001), 
respectively (Fig. 2a,b, Table 2). 

Differences between the arable controls and bioenergy 
crop were smaller than those seen in the grassland tran- 
sitions, with greater variation between sites, and not sig- 
nificant (P = 0.071 and P=0.846 for SRC and 
Miscanthus transitions, respectively) (Fig. 2). The non- 
significant differences in mean soil C stocks are being 
16.27 + 7.18 Mg C ha^! higher under SRC, and 
2.26 + 8.18 Mg C ha^! lower under Miscanthus planta- 
tions compared to arable controls. 

Within the SRC data, the grassland control at site 17 
had exceptionally high soil C compared to its paired 
bioenergy crop (Fig. 2a). This transition unit was 
located at a site with highly complex underlying geol- 
ogy and variable soil types. Removing this transition 
from the analysis of soil C stock reduced the difference 


between the SRC and the grassland control. The mean 
soil C stock under the SRC, however, was still signifi- 
cantly lower (—23.34 + 8.37 Mg C ha‘ ') than the grass- 
land controls (P = 0.047). 

Differences in soil C stocks between the bioenergy 
crops and the controls were reflected in the annual rates 
of change (Mg C ha! yr )) in the surface soil (0- 
30 cm) with effects of both the original land use 
(F157 = 11.99, P = 0.001) and also bioenergy crop type 
(Fi57; = 6.59, P = 0.014) but there was no interaction 
between these factors (F; 5; = 0.326, P = 0.571) (Table 3). 
Rates of change in the transitions from grassland, as 
would be expected by the differences in soil C stock, 
were consistently negative and significantly lower than 
observed in the arable transitions. Unlike the differences 
in soil C stock, annual rates of change also allowed the 
comparison of the two bioenergy crops and showed that 
the rates of change for the SRC transitions were more 
positive that those for the Miscanthus transitions 
(Table 3). 


Soil C stocks 0—1 m 


Over 0-1 m, soil C stocks (Mg C ha!) and annual rates 
of change followed a similar pattern to those seen in the 
surface soils (Fig. 1c,d, Tables 2 and 3). Unlike the sur- 
face soil, however, differences in soil C stocks between 
the controls and bioenergy crops were not significant in 
either the SRC GZ (3) = 1.93, P = 0.3813, R2 = 0.92) or 
Miscanthus transitions D? (3)-22.10, P = 0.350, 
R — 0.90)] (Table 2). Annual rates of change were not 
significantly different between the bioenergy crops 
(F134 = 0.015, P = 0.902), nor was there any impact of 
the original land use (Fi34,- 2.432, P = 0.128) or an 
interaction between these factors (F134 = 1.166, 
P = 0.287) (Table 3). 

Over a shallower depth of 0-50 cm, there were differ- 
ences in soil C stocks in the SRC transitions (7? 
(3) = 7.16, P = 0.028, R2 = 0.91) but not the Miscanthus 


Table 2 Mean litter and soil C stocks (Mg C ha ') and standard error for the bioenergy crops (SRC and Miscanthus) and controls 


C stock (Mg Cha" 


Land use Litter 0-30 cm 0-50 cm 0-100 cm 

SRC 0.97 + 038? 70.31 + 6.57? 91.16 + 8.98? 116.91 + 11.65? 
Arable 0.38 + 0.19^ 54.04 + 8.18? 76.41 + 12.14? 107.22 + 16.14? 
Grassland 0.21 + 0.19» 103.87 + 9.5^ 129.03 + 12.6* 147.19 + 16.56? 
Miscanthus 2.09 + .0.24? 74.31 + 7.84? 108.77 + 7.70? 124.31 + 11.39? 
Arable 0.78 + 0.33? 76.57 + 8.53* 94.41 + 9.51? 120.98 + 12.69? 
Grassland 0.06 + 0.36^ 101.14 + 8.895 12347 + 10.14? 140.49 + 13.75? 


0-50 cm ESM and 0-100 cm ESM refer to soil C stock based on reference soil mass for these depths of 6 and 13 Gg ha 1 Same litter 
indicates nonsignificant difference > P 0.05; * indicates that there was a near-significant difference (P = 0.063) between the grassland 


and the SRC. Test conducted on Miscanthus and SRC transitions separately and within each depth division. 
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Table 3 Annual rates of change in soil C stocks for 0-30 cm, 0 
—50 cm and 0-1 m soil cores based on ESM. Annual rates of 
change are estimated by dividing change in mean soil C com- 
pared to control by the years since transition. n = 15 and 3 for 
the 30-cm cores and 1-m cores, respectively 


Rate of change Mg C ha“! yr! (SE) 


Land-use 

Change 0-30 cm 0-50 cm 0-1 m 

SRC vs. 1.54 + 0.70 1.93 + 1.37 1.26 + 1.41 
Arable 

SRC vs. 1.69 + 0.81 2.98 + 1.61 2.74 + 1.65 
Grassland 

Miscanthus —0.93 + 0.74 0.18 + 1.37 0.05 + 1.41 
vs. Arable 

Miscanthus 3.17 + 0.81 2.11 + 1.52 0.69 + 1.65 


vs. Grassland 


transitions (y* (3) = 4.34, P = 0.114, R2 = 0.85) (Table 2). 
The significant difference in the SRC transitions was, 
however, related to differences in soil C stocks between 
the grassland and the arable control (P = 0.008), 
although there was also nonsignificant trend for lower 
soil C stocks within the grassland controls compared to 
the SRC (P = 0.063). 

Rates of change reflected the absence of a significant 
difference in soil C stock, which were similar in both 
bioenergy crops (F135 = 0.188, P = 0.667). There was 
no interaction between the current land use and the 
control land use (F135 = 0.761, P = 0.388) but rates of 
change were lower in the grassland compared to ara- 
ble transitions (F5 = 5.952, P = 0.019; Table 3), high- 
lighting a difference that was less clear with soil C 
stock. 


Driving factors determining changes in soil C 


Based on the model selection, soil C stocks of the con- 
trol field and the current land use (SRC, Miscanthus) 
were tested for their effect on the percentage change in 
soil C stocks resulting from the transition the bioenergy 
crops (Tables 53, S4). Soil C stocks was found to be neg- 
atively related to the percentage difference in soil C in 
bioenergy fields (? (1) = 870, P = 0.003, R2 = 0.52). 
There was no interaction between current land-use type 
(SRC, Miscanthus) and soil C stock (Gi (1) = 2.138, 
P = 0.144, R = 0.51), suggesting a similar relationship 
in both SRC and Miscanthus, but a near-significant effect 
of land-use type was observed Ca (1) = 3.216, P = 0.073, 
R2 — 0.22) likely resulting from the different intercepts 
of the linear relationships in the two bioenergy crops 
(Fig. 3 a & b). Examination of the residuals highlighted 
that three transitions (Miscanthus transitions 24 and 25, 


and SRC willow transition 17) had a large influence on 
the results. Removal of these transitions influenced the 
slope of the linear relationships (Fig. 3 c & d), but did 
not change the overall significance of any of the factors. 

Time since bioenergy establishment and the clay con- 
tent of the bioenergy crop were the third most impor- 
tant factors influencing the change in soil C stock based 
on the marginal R? and AIC scores, respectively 
(Table S3). However, there was no clear relationship 
between the percentage change in soil C stock and 
either time since bioenergy establishment or clay con- 
tent (Fig. 4). 


Litter C stocks 


Litter C stocks were different between the land uses in 
both the Miscanthus (7 (2) = 25.42., P = 0.001, R = 0.84) 
and the SRC plantations (7 (2) = 43.68, P < 0.001, 
R — 0.69), with post hoc testing showing that litter stock 
was higher in the bioenergy crops than in either the ara- 
ble or grassland controls (Table 2). The addition of these 
relatively small litter C stocks to the surface soil C stocks 
(0-30 cm, Table 2) has little effect on the impact of the 
bioenergy crops on C stocks. C stocks remain lower in 
the bioenergy transitions than in grassland controls and 
are not significantly different to the arable controls. 


Discussion 


Soil C stocks in arable transitions 


In this study, annual rates of change in the surface soil 
were more positive for the arable transitions than for 
the grassland transitions. Although, as soil C stocks in 
the SRC and Miscanthus plantations were not signifi- 
cantly different to the arable controls, the difference in 
the rates of changes is most likely related to the nega- 
tive impacts on soil C stocks of transition from grass- 
land, rather than any positive impacts of arable. This 
absence of a positive impact is contrary to a number of 
studies which have reported increases in topsoil C 
stocks following transitions from arable land uses to 
these bioenergy crops (Jug et al., 1999; Dondini et al., 
2009; Schmitt et al., 2010; Felten & Emmerling, 2012). 
These studies used a fixed depth method (FD) to calcu- 
late soil C stock which, unlike the ESM used in this 
study, makes no adjustment for changes in bulk density 
(BD) (Bárcena et al., 2014). Applying FD methods to our 
data leads to a similar result to these studies with signif- 
icantly lower surface soil C stock in the arable controls 
(Tables S5 and S6). The use of a FD method appears to 
inflate the differences between the arable control and 
the bioenergy crops, something that has been noted in a 
similar land-use change study (Bárcena et al., 2014). The 
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Fig. 3 Relationship between 0-30 cm soil C stocks in control crops and percentage differences for control in soil C resulting from 
land-use change for: SRC transitions (a), Miscanthus transitions (b), SRC transition without site 17 (c), Miscanthus transition without 
transitions 24 and 25 (d). Red markers indicate arable transition green grassland transition. The line shows linear regression of change 
in soil C stock with C stocks of the control fields; shaded area shows 95% of confidence interval, R? gives values for individual regres- 


sion lines. 


use of an ESM method is not widespread in bioenergy 
studies, and in the case of arable transitions, the only 
comparable study is that by Walter et al. (2014). Using 
an ESM method and a paired-site approach to assess 
impacts of arable to SRC transitions, Walter et al. (2014) 
also reported consistent changes in surface (0—30 cm) 
soil C stock. 

Below the plough layer, BD is more consistent 
between land uses, and differences in C stock estima- 
tion due to method are less apparent. This is possibly 
reflected by the studies that have assessed soil C stock 
below 30 cm and reported no significant changes in 
transitions to either SRC (Coleman et al., 2004; Lockwell 
et al., 2012; Bonin & Lal, 2014; Walter et al., 2014) or Mis- 
canthus (Felten & Emmerling, 2012). 

The age of the plantations studied may also have an 
impact on the soil C stock change. Hansen et al. (2004) 
reported higher soil C stocks under Miscanthus planta- 
tion compared to arable controls but only under the 
older of two plantations sampled (9 and 16 years old). 
A study of SRC by Dimitriou et al. (2012) also reported 
an increase in soil C stock compared to arable controls, 
but only one of the 14 sites sampled was under 15 years 
old. In addition, many were not in optimum condition 


leading the authors to suggest that some of the increase 
in soil C concentration could be related to C inputs from 
decaying stools and roots. 

Within this study, the difference in the mean age of 
the bioenergy crops may also explain the differences in 
the rate of change between the SRC and the Miscanthus 
transitions. The mean annual rate of change in the sur- 
face soil for the SRC to arable transitions 
(143 + 0.71 Mg C ha! yr D, with a mean age of 
8.5 years, was within the upper range of reported val- 
ues from 0.38 to 1.59 Mg C ha! yr ! (Kahle et al., 2010, 
2013; Chimento et al., 2014). In contrast, the annual rate 
of change for transitions to Miscanthus from arable 
(—0.93 + 0.74 Mg C ha^! yr ), with a mean age of 
6.4 years, was more negative than the mean reported 
values for topsoil changes of 0.28-2.24 Mg C ha^! yr! 
(Clifton-Brown et al., 2007; Dondini et al., 2009; Zimmer- 
man et al., 2012; Chimento et al., 2014). This possibly 
reflects the mature plantations in some of these studies 
(16 years and 14 years in Clifton-Brown et al., 2007 and 
Dondini et al., 2009; respectively) compared to this 
study. It is also clear that impacts on soil C vary greatly 
between sites, even within individual studies. For exam- 
ple, although mean rates of change in the study by Zim- 
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Fig. 4 Relationship between time since establishment (a) and bioenergy clay content, (b) and the percentage change in 0-30 cm soil 
C. Red markers indicate arable transition green grassland transition, diamond indicates SRC transitions and squares Miscanthus transi- 


tions. Error bars show pooled SE. 


merman et al. (2012) were 1.72 Mg C ha! yr" from 
arable, the rates of change across sites within this study 
ranged from —6.85 to 7.7 Mg C ha ! yi Site-specific 
factors clearly influence the impacts on soil C stocks, as 
reflected in the between-site variability observed within 
this study and reported in other multi-site studies 
(Coleman et al., 2004; Dimitriou et al., 2012; Don et al., 
2012; Walter et al., 2014). 

One possible additional source of variability between 
sites could be related to the willow clones selected. 
Nearly all the sites visited were planted by a single con- 
tractor whose records do not contain details of the 
clones planted at each site (F. Walters, Coppice 
Resources Ltd, Retford, pers.com.) but only that the 
mixed willow will contain 4—5 different clones. The lack 
of detailed information coupled with the practice of 
mixing clones throughout a single plantation (e.g. clones 
are not planted in uniform strips) for pest control pur- 
poses means that it is not possible within this study to 
examine differences between the influence of individual 
clones. However, any differences in soil C stock result- 
ing from different clones are likely to be smaller than 
the impact resulting from the LUC from arable or grass- 
land land uses. 


Soil C stocks in grassland transitions 


In contrast to the findings for arable soils, the lower soil 
C stocks in the topsoil (0-30 cm) and the negative rates 
of change of the SRC and Miscanthus plantations com- 


pared to the grassland controls reflect findings in other 
studies (Don et al., 2012; Rytter, 2012; Zimmerman et al., 
2012). The mean annual rate of change in the transition 
to SRC (—1.69 + 0.82 Mg C ha! yr !, 0-30 cm) com- 
pares well, once again, with the values reported in a 
study of a 9-year-old SRC willow plantation by Lock- 
well et al. (2012) of —222 and -111 Mg C ha! yr! 
over 0-20 cm and 0-40 cm depths, respectively. The 
mean rate of change for transitions to Miscanthus from 
grassland, however, was again more negative 
(—3.17 + 0.81 Mg C ha! yi D than those reported 
from —1.66 and 0.83 Mg C ha! yr! by Zimmerman 
et al. (2012) and Zatta et al. (2014). 

Over the greater depth of 1 m, the magnitude of 
differences in soil C stocks observed was similar to 
those seen in the surface soil, especially for the SRC 
transitions (—33.55 + 7.52 Mg C ha! and —30.28 + 
10.96 Mg C ha^! for 0-30 cm and 0-1 m, respectively) 
but differences were no longer significant. Fewer 1-m 


cores were taken compared to the 30-cm cores, resulting 
in reduced statistical power to detect impacts at greater 
depth. Walter et al. (2014) and Lockwell et al. (2012), 
however, reported similar findings in transitions from 
grassland to SRC, concluding that soil C losses in the 
surface soil were offset by increases lower in the soil 
profile, resulting in no significant changes in soil C 
stocks overall. Miscanthus shares the tendency of SRC to 
be deep rooting and Miscanthus-derived C inputs have 
been detected at depths of up 1.5 m (Felten & Emmer- 
ling, 2012), and thus, there is a mechanism by which 
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both crops could alter soil C stocks at depth. In this 
study, mean difference in soil C stocks between both 
the SRC and Miscanthus and their grassland controls 
was less negative over 0-1 m than over 0-30 cm. Differ- 
ences in sampling intensity between 0-30 cm and 0-1 m 
cores mean that it is not possible to directly attribute 
any redistribution of soil C within the soil profile. 

Alternatively to a redistribution of soil C stocks, it is 
possible that changes in soil C stock were limited to the 
surface soil and that difficulties in detecting changes in 
soil C stock 0-1 m are instead due to the dilution of the 
impacts in the surface soil when including soil C stock 
at greater depths. This would agree with studies which 
report slower turnover times in the subsoil, with 
reported mean C resident times in soil layers below 
20 cm of 2000-10 000 years (Fontaine et al., 2007). Sam- 
pling subsoil is, however, still extremely valuable as 
although C stocks at depth may be characterised by 
long residence times, they have also been found to be 
susceptible to priming resulting from labile C inputs 
such as root exudates (Fontaine et al., 2007; De Graaff 
et al., 2014). Deep soil coring therefore provides a mech- 
anism to detect both increase in soil C and any losses 
due to C priming. 

Regardless, if losses in the surface soil are replaced 
with gains at depth or just diluted, any step taken to 
reduce surface soil C loss would be beneficial. Grass- 
land soil C stocks have been shown to be negatively 
affected by tillage (Poeplau & Don, 2014). Thus, it has 
been suggested that the intensive cultivation under- 
taken prior to the bioenergy crop establishment may 
account for a substantial proportion of soil C losses 
observed (Don et al., 2012; Walter et al., 2014). A move 
to new, less intensive establishment methods may pro- 
vide one option to reduce impacts on soil C stocks. 
However, it is unclear what role other factors, such as 
changes in the quality or quantity of inputs to soil, may 
play in addition to the effects of cultivation. For exam- 
ple, Poeplau & Don (2014) reported that transitions 
from grassland to forestry resulted not only in changes 
in soil C stocks but also a shift in soil C from stable to 
labile pools. 


Factors influencing changes in soil C stock 


Explaining variations in soil C stock changes within this 
study was explored through assessment of relationships 
between changes in soil C stock and selected factors. A 
negative relationship was found between changes in soil 
C stock and the soil C stock of the control field, suggest- 
ing that establishment of bioenergy crops on sites with 
low initial soil C provided the best opportunity to 
derive positive impacts on soil C stocks. Such a negative 
relationship was predicted for SRC poplar plantations 


in modelling work by Garten et al. (2011) and generally 
agrees with the conclusions of Don et al. (2012) and 
Walter et al. (2014) that conversion of arable lands, 
which generally have low soil C stock, is preferable to 
conversion of grassland for bioenergy plantations. It is 
difficult to separate the impacts of original soil C stocks 
and original land use because they are highly correlated 
(e.g. higher soil C stocks are generally associated with 
grassland sites). As land use also affects soil C stability 
and turnover, as well as soil C stocks (Poeplau & Don, 
2014), impacts of land-use transitions could be influ- 
enced by both the stability of the soil C and the total soil 
C stocks. 

The relationships between control soil C and changes 
in soil C stock following bioenergy crop establishment 
are relatively weak, especially for Miscanthus. Transi- 
tions to SRC and Miscanthus have R? values of 0.30 and 
0.01, respectively, which indicate considerable unex- 
plained variability related to the impacts on soil C 
stocks at individual sites. Part of this unexplained vari- 
ability may reflect the challenge of finding paired sites 
with no pre-existing differences in soil C stocks between 
the two land uses before conversion. In many cases, the 
bioenergy crops and paired sites were adjacent but the 
soil texture analysis does suggest, even for the more 
reliable sand and silt data, that in a few of the sites 
there may be some underlying differences between 
some of the transition pairs. In addition, whilst finding 
sites with generally similar land-use histories was rela- 
tively straightforward, the normal crop rotation prac- 
tices (rotations wheat, barley, beans, etc.) and the 
variable nature of farming (fertiliser inputs, harvest 
times, etc.) combined with the limited nature of long- 
term data held by land owners meant that some vari- 
ability between the bioenergy crop and the paired con- 
trol was inevitable. A better understanding of between- 
site variability is also clearly needed. For example, in 
this study, the rates of change in soil C stock range from 
—3.75 to 0.58 Mg C ha’! yr! for grassland to SRC tran- 
sitions, and from —7.44 to 2.53 Mg C ha! yr! for 
grassland to Miscanthus transitions. 

It is worth noting that whilst underlying differences 
between the paired sites could have influenced the anal- 
ysis of the potential factors driving soil C stock change, 
and the rates of change, where the percentage change 
was calculated at the transitions level, in the assessment 
of soil C stocks individual core data rather than transi- 
tion level mean were used. When using this core data, 
the mixed model is less sensitive to variation between 
the bioenergy crop and the control. 

No relationship was found between bioenergy plana- 
tion age or clay content and changes in soil C stock. In 
the case of clay content difficulties with both the analy- 
sis method and a limited range of clay content across 
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the sites (3.50-12.56% Table 1 Fig. S1) may have 
reduced our ability to detect a relationship. However, 
the absence of any relationship between soil texture and 
changes in soil C stocks has also been reported for SRC 
(Walter et al., 2014) and Miscanthus (Poeplau & Don, 
2014). Clay content tends to be positively associated 
with soil C stock (Stockmann et al., 2013) and the 
absorption of C compounds to clay minerals, together 
with occlusion into clay aggregates, has been shown to 
stabilise soil organic matter (Dungait et al., 2012; Stock- 
mann et al., 2013). Therefore, there could be an expecta- 
tion that higher clay content would protect soil C 
during LUC, and aid its accumulation post LUC (Laga- 
niere et al., 2010). One possible reason why this is not 
seen could be that the current practice of intensively till- 
ing prior to bioenergy crop planting could reduce the 
protection afforded by occlusion into clay aggregates 
(Stockmann et al., 2013). 

A relationship between changes in soil C stock and 
time since bioenergy crop establishment was also 
absent, something which has been reported in a num- 
ber of other multi-site studies (Don et al., 2012; Walter 
et al., 2014). This is despite general agreement across a 
wide range of land-use transitions that time since LUC 
is an important factor in determining soil C stocks 
(Bárcena et al., 2014; Poeplau & Don, 2014; Walter et al., 
2014). Bárcena et al. (2014) suggested that the time 
taken for soil C stocks to recover from any initial soil 
C loss following land-use transitions, and to reach a 
new equilibrium, may vary between sites. Thus, any 
assessment made between sites that are yet to near a 
new equilibrium will lead to highly variable results 
(Bárcena et al., 2014). In case of transitions from arable 
to forestry, Bárcena et al. (2014) found that increases in 
soil C were only detectable in a chronosequence of 
independent sites after 30 years. The time required for 
soil C recovery in SRC and Miscanthus plantations is, 
as yet, unknown. Walter et al. (2014) did select older 
plantations (15-35 years) in their study of 21 SRC plan- 
tations, but were still unable to detect any relationship 
between plantation age and impacts on soil C stock. 
Therefore, it may be that the time period required to 
detect an effect of age on soil C under bioenergy crops 
will exceed the expected 25-30 year life span of these 
plantations. 

The time taken to reach a new soil C equilibrium has 
potential to impact on the ^payback time' required for 
any decreases in soil C stocks within the soil to be 
replaced (Mello et al., 2014). In contrast to the transi- 
tions from arable, where changes in soil C stock were 
not significant, there is not a soil C debt to be paid. A 
soil C debt was detected in the surface soil, at least in 
grassland transitions. To replace this debt through 
increases in the soil C stock, the bioenergy crops must 


in theory reach a new soil C equilibrium that is equal to 
or greater that than of the grassland. The time it takes 
to reach this new equilibrium is also critical because, if 
it takes longer than the lifetime of the bioenergy crop, it 
may not be possible to repay the soil C debt through 
changes in soil C stock alone (Bárcena et al., 2014; Mello 
et al., 2014). Although it must be recognised that over 
greater depths this and other studies have found no sig- 
nificant negative impact of planting on grassland (Wal- 
ter et al., 2014). Although requiring a detailed life-cycle 
assessment to confirm, the C saving attributed to using 
biomass to offset fossil fuel use may be greater than any 
soil C loss as has been found to be the case for sugar- 
cane planted on pasture in Brazil (Mello et al., 2014). 

It is possible that the difficulties in detecting a clear 
chronosequence may also result from different sites 
having different linear relationships between age and 
soil C and/or more complex nonlinear relationships. In 
addition, the C stock within the control field may not be 
in equilibrium, and for this reason, it is best to view 
controls as counterfactuals rather than a time zero. 
Long-term studies utilising both repeated sampling and 
the use of counterfactual paired sites, soil fractionation 
(Poeplau & Don, 2014) and process base modelling 
(Dondini et al., 2015) are all methods which could help 
to provide a better understanding of the time it will take 
to reach a new equilibrium, and allow the comparison 
to other land-use options. The data collected in this 
study are highly suited for process models, which can 
be used to understand key drivers of soil C change, and 
such models can be used to predict impacts of future 
climate scenarios (Dondini et al., 2015). 

We conclude that where choices exist, the selection of 
arable land for bioenergy transitions to SRC and Miscant- 
hus is likely to be more positive for soil C stocks than 
conversion from grassland, at least for soil C stocks 
within the surface soil. Whilst changes in soil C stocks at 
0-1 m were not significant in any of the transitions types, 
the direction of changes mirrored those in the surface 
soil. Questions still remain as to why transitions from 
grassland can lead to negative changes in soil C, and 
work on soil C stability, especially during bioenergy crop 
establishment, would both address this question and 
potentially provide insight into management solutions 
that would maximise the soil C sequestration potential of 
these crops. Whilst these conclusions are valid for soil C, 
the findings also need to be considered in the wider con- 
text of other ecosystem services such as productivity, 
greenhouse gas regulation and water quality. 
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